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Abstract 

We examine the leading order semiclassical gradient corrections to the non-interacting kinetic energy 
density functional of a two dimensional Fermi gas by applying the extended Thomas-Fermi theory at finite 
temperature. We find a non-zero von Weizsacker-like gradient correction, which in the high-temperature 
limit, goes over to the functional form (/i^/24m)(Vp)^//9. Our work provides a theoretical justification for 
the inclusion of gradient corrections in applications of density-functional theory to inhomogeneous two- 
dimensional Fermi systems at any finite temperature. 

PACS numbers: 71.15.Mb,03.65.Sq,05.30.Fk,31.15.Bs 
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I. INTRODUCTION 



In 1966, Power et al- performed transport measurements on a Si-metal-oxide-semiconductor 
structure in which a degenerate gas of electrons was electrostatically induced. Their work demon- 
strated for the first time that the density of states in the n-type electron inversion layer had the 
expected behaviour for a two-dimensional electron gas (2DEG). Since Fowler's seminal work, 
the exploitation of the electronic properties of III-V semiconductors has led to the realization of 
high quality, high mobility 2DEG's at the interface of epitaxially grown III-V structures, such 
as GaAs/AlGaAs heterostructures.- Through electrostatic, and/or etching techniques, the 2DEG 
found in the III-V semiconductor interface, can be manipulated to create experimental realizations 
of low dimensional electron systems such as quantum wires, quantum dots and quantum anti-dots. 

By far, the workhorse for a theoretical understanding of the bulk electronic properties of such 
low-dimensional electronic systems is the zero-temperature (T = 0) density-functional theory 
(DFT) of Hohenberg, Kohn and Sham (HKS).^ '^ The key element in the HKS approach is the defi- 
nition of the kinetic energy (KE), corresponding to a system of noninteracting fermions moving 
in some effective one-body potential. The HKS scheme treats the KE exactly at the independent 
particle level, and so the development of explicit, orbital-free functional for T[p\, the KE den- 
sity functional, is an important objective. Ideally the appropriate functional should yield both the 
correct energy and the correct density profile. 

To this end, the simplest approach for the construction of the KE density functional, T[p], is the 
local-density approximation (LDA), sometimes referred to as the Thomas-Fermi (TF) approxima- 
tion.'^ In this approximation, the known form for KE density functional of the uniform electron 
gas, is also used locally for the KE density functional of the inhomogeneous system. One would 
then expect the LDA to be applicable only in cases where the density in the system is a slowly vary- 
ing function of position. In fact, this is not the case, and even in highly inhomogeneous systems, 
the LDA is found to work reasonably well.^'^ 

Although the LDA for the KE density leads to reasonable results for the energy, the calculated 
density profile in a self-consistent DFT scheme does not exhibit the desired quantum mechan- 
ical tunnelling into the classically forbidden region. To overcome this issue, the so-called von 
Weizsacker (vW) gradient correction,^ oc (Vp)^/p, is added to the KE functional. In 3D, the vW 
gradient correction can be rigorously justified within the extended TF (ETF) theory, originally de- 
veloped in the context of nuclear physics.^^^ The inclusion of the vW term leads to smooth and 
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continuous densities, while improving the quality of the KE functional by taking into account the 
inhomogeneity of the system. 

An application of the ETF theory to 2D systems, however, leads to the conclusion that there are 
no gradient corrections to the 2D KE density functional.— — This, of course, makes no physical 
sense since the LDA cannot be variationally exact for an inhomogeneous system. Thus in DFT 
applications to systems derived from the inhomogeneous 2DEG discussed above, a phenomeno- 
logical approach must be taken in which a vW-like gradient correction is put in "by hand".— 
Although the vW-like correction term is entirely ad hoc for a 2D system, its use has been justified 
by (i) the KE reduces to the TF limit for slowly varying densities, and (ii) it allows one to represent 
strongly inhomogeneous densities in a quantum mechanically reasonable way. 

To date, there has been no formal justification for the inclusion of a vW-like term for 2D sys- 
tems at zero-temperature. In this paper, we establish the existence of a vW-like gradient correction 
to the 2D KE density functional at finite temperatures. Our approach parallels the earlier work of 
Brack/ ^ in which the ETF theory was developed in the context of "hot" nuclear matter (ETFT). 
Given the recent work of Eschrig,— which aims at providing a rigorous foundation for DFT at 
finite-temperature, the results presented in this paper are immediately relevant to future applica- 
tions of T > DFT in low dimensional electronic systems. 

The rest of our paper is organized as follows. In the next section we provide a brief review 
of the general ETFT approach, followed by an explicit calculation of the T > second order 
gradient correction to the 2D KE functional. In Sec. Ill we numerically investigate the quality 
of the gradient corrected functional by comparing it to known, exact results, for an isotropic 2D 
harmonic oscillator at finite temperature. The paper concludes in Sec. IV with a brief summary 
and suggestions for future investigations. 

II. EXTENDED THOMAS-FERMI THEORY AT FINITE-TEMPERATURE 

In this section, we provide a brief review of the ETFT approach. The interested reader should 
refer to the Ref. Ill ill for a detailed discussion of the ETFT and the Wigner-Kirkwood semiclassical 
expansion. 
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A. Semiclassical spatial density 



At the heart of the ETFT approach is the Wigner-Kirkwood (WK) semiclassical expansion of 
the zero-temperature, diagonal Bloch density matrix (BDM), which in 2D is given hy^ 



a(r;/3) 



-/3V(r) 



12m 



VV 



(3 



+ 



(1) 



where V{r) is a local one-body potential, and for our purposes, we have only shown terms up to 
relative order fi^ and A = (27r/i^/3/m)^/^. Note that (3 is to be viewed as a complex parameter here, 
and not the inverse temperature 1 / (/^sT). In order to incorporate finite temperatures into the WK 
semiclassical theory, the finite-temperature BDM is defined by^ 



CT(r;/3)^Co 



(2) 



sm{TT/3kBT) 

The finite-temperature spatial density is then obtained from an (all two-sided) inverse Laplace 
transform (ILT) of the finite-temperature BDM, viz., 



(3) 



where the factor of two in Eq. (3) accounts for the spin degeneracy after the spin trace has been 
taken, and /i has the physical significance of the chemical potential. It should be noted that if the 
exact Ct{t; (3) is known, then Eq. (3) will yield the exact, quantum mechanical, finite-temperature 
spatial density. Of course, here, we are using a semiclassical expansion for CT(r;/3), so the 
resulting p(r; T) will be the semiclassical spatial density. 
In what follows, we will make use of the following ILTs: 
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where Li2(-) is the polylog function,— and z = exp[(/i — V)/kBT]. Using Eq. (3), along with 
Eqs. (4-6), readily leads to the following second-order expression for the finite-temperature spatial 
density: 



P(r;T) 
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The ETFT density in Eq. (7) is well-defined throughout all space, with the last two terms, denoted 
by p'^^^(r; T), being relative order h? greater than the first term, p^^\r] T). Note that in Eq. (7), for 
V(r) < /i, 2; — 7- oo exponentially as T — 0, so that the T -> limit is well-defined only within the 
classical region; the non-analytic behaviour of the zero-temperature ETF densities at the turning 
point, /i = V(r), is well-known. What we find here is that the singular behaviour of the T = 
densities cannot be avoided by first formulating the semiclassical theory at finite-temperature, and 
then performing the T — )■ limit.'- We have, however, confirmed that the T — )• limit of Eq. (7) 
with V{y) < ^ correctly reduces to the known (albeit problematic) T = result.— 



B. Semiclassical kinetic energy density 

The KE density may be obtained from knowledge of the finite-temperature first-order density 
matrix (EDM). To this end, it is useful to introduce the centre-of-mass, q = (r + r') / 2, and relative 
coordinates, s = r — r', so that we may write three variants of the KE density I'^Si^i 

(\ 



r(r; T) = -— (^-V^ + ) p(q, s; T) |,=o , (8) 
ri(r; T) = ^ fiv^ - ) p(q, s; T) |,=o , (9) 



2m V4 
r(r;T) + ri(r;r) 

Ur;T) = ^ = -— V,p(q,s;r)|s=o . (10) 

Again, if the exact finite-temperature expression for the EDM is known, then Eqs. (8-10) will 
yield the exact, quantum mechanical finite-temperature KE density. 

While all three of the above expressions for the KE density integrate to the exact same kinetic 
energy, 7i(r) is strictly positive definite, and is therefore sometimes preferred in applications of 
density functional theory. It has already been shown long ago that T(r; T) and 7i(r; T) generally 
have oscillations exactly opposite in phase, so that their mean, ^(r; T), is a smooth function. In 
this paper, we focus on T(r; T) in order to make contact with earlier theoretical work done at 
zero-temperature, where the exact T(r; T = 0) was investigated (see also Eq. (26) in Sec. III).— 

The kinetic energy density in Eq. (8) may also be expressed in terms of only local quantities, 
viz.^ 

r(r; T) = -p(r; T)\/(r) + ^(r; T) + Ta(r; T) , (11) 

where 



i?(r; T) = pp(r; T) - 
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is the free energy density, and 



a{r-T) = —C^ 



dT' 



Cn 



(13) 



/32 

is the entropy density. The semiclassical approximation to T(r; T) may easily be determined to 
second-order by employing the semiclassical approximation to C^ij] (3), as in Sec. IIA for the 
spatial density. A straightforward calculation results in the following expression for the finite- 
temperature, semiclassical KE density 



nh? ' ' l2-KkBT{z + 1)2 1271 z + l 
= r^'>\r-T) + T^^\r-T) , (14) 

where, following Eq. (7), the last two terms in Eq. (14) are denoted collectively by T''^^(r; T). As 
in Eq. (7), the T — )• limit of Eq. (14) is well-defined only within the classical region. Finally, it 
follows immediately from Eqs. (8) and (9) that 

ri(r;T) = r(r;r) + |?-VV. (15) 
C. Second-order kinetic energy density functional 



For the special case of 2D, the elimination of z and V in T(r; T) above, in favour of the spatial 
density, p, is quite straightforward.-^2i2^'2^ To begin we define V = so that z = exp(V^). 

Thus, the density, Eq. (7), is a function of V, VV and V^V. Calculating from Eq. (7) Vp and 
V^p and consistently neglecting higher than second derivatives of the potential, we have: p = 
p{V,VV,V'^V), Vp = Vp(- ■ and VV = VV(- ■ ■)' which can be solved for V,VV and 
V^V^; inserting this into Eq. (14) yields the finite-temperature kinetic energy density functional up 
to Oih^), viz., 

TetftH = -ATfcBTLi2(l-e^/^-)--^VV-T^/i(r;T)VV 



12m 12m" 



where 



24m p 

TtftW+Te^ftW, (16) 
/2(r;T) = e^/^-[/i(r;T)]2, (18) 



and At = mkBT / (Tih?). The first temi in Eq. (16) is the finite-temperature 2D TFT KE density 
functional, 7tft[p], while the other three terms represent the 0{h^) gradient corrections. As 
advertised, the last term in Eq. (16) has the vW form ~ ( Vp)^/ p. 

An explicit expression for the 7i[p] KE functional may also be given by making use of Eq. (15), 
viz., 

ri,ETFT[p] = -ATkBTLi^il - eP/^n + I^VV - -^fi{r; T)V'p 

bm 11m 

z4m p 

Recall that 7etft[p] and 7i,etft[p] both integrate to the same total kinetic energy for finite sys- 
tems since the Laplacian term is the divergence of a vector field which vanishes at infinity, and by 
Gauss's theorem will not contribute the kinetic energy. 

To investigate the low temperature behaviour of Eqs. (16) and (19) we use, 

lim U,{-e^l^-) = -\(^)\ (20) 

along with the fact that /i(r; T) 0, and /2(r; T) as T 0, to write 

TetftIp] ^ ^ (^vrp2 _ IvV^ , (21) 

and 

Ti.etftM ^ 1^ (^V + ^VV^ . (22) 

Equations (21) and (22) agree with the known results for the T = 2D KE functionals T[p] and 
7i[p], respectively.^^ As mentioned above, for physical densities, integration over V^p vanishes. 
Therefore, in any practical implementation of self-consistent T = DFT, the Laplacian term may 
be ignored, and we may write 

TetftM = TI.etftIp] = — (ttp^) , (23) 

as T — ;> 0. We see that there is no vW-like gradient correction at T = 0, leading to the incorrect 
conclusion that for an inhomogeneous 2D Fermi gas, the TF KE functional (at least to second- 
order) is exact.^^ We would like to stress again that the non-uniqueness of the KE density does 
not alter the result that there are no vW-like gradient corrections at T = 0; the only differences 
between the T = semiclassical KE densities obtained from Eqs. (8-10) are Laplacian terms 
which, as we have already stated, are of no consequence since they vanish upon integration for 
physical (i.e., finite) systems. 
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It is readily found that as T ^ oo, /i(r; T) 1 and /2(r; T) 1, so that up to C>(/i^), the 2D 
KE functionals go over to 

Tetft [p] ^ Tb [p] = A:bTp - ^ V V + , (24) 

6m 24m p 

and 

ri,ETFT[p] ^ T,,b[p\ = ksTp + -^VV + 7^^^ . (25) 

12m 24m p 

Therefore, at high-temperature, the second order gradient corrections take a functional form anal- 
ogous to what is found in 3D ETF, although the numerical pre-factors are different. It is also 
interesting to note that the high-temperature limit of the TFT term is linear in the density, ksTp, 
in contrast to the quadratic dependence, 7rp^/2, exhibited in the zero-temperature limit. We have 
also checked that the high-temperature limit of the 2D KE functional may be obtained by perform- 
ing an analogous calculation assuming a Boltzmann, rather than a Fermi, gas. 

ni. COMPARISON WITH EXACT RESULTS 

In a previous study. Brack and van ZyP^ examined the T = 2D TF KE functional by com- 
paring its global {i.e., integrated) and local {i.e., spatially dependent) properties with the known 
analytical expressions for the 2D harmonic oscillator (HO) potential. They found the remarkable 
result that at T = 0, the 2D TF KE functional {without gradient corrections), when using the 
exact spatial density of the 2D HO, gives the exact quantum mechanical kinetic energy. This result 
is highly non-trivial because the TF functional is simply the LDA to the true KE, and therefore, 
cannot be variationally exact. More surprising, however, is how well the local behaviour of the 



exact KE density is reproduced by the TF approximation, as illustrated in Fig. 3 of Ref. 112211 . The 
purpose of this section is to perform an analogous calculation for the finite-temperature KE den- 
sity functionals presented in this paper. In our numerical calculations, we have scaled all energies 



and lengths by fku, and £osc = \/h/muj, respectively. We have also restricted our attention to 
relatively small particle numbers, N, since it can be shown rigorously that in the large- limit, 
the TF approximation becomes exact.— 

The exact finite-temperature kinetic energy density is given by (see also Eq. (8)) 

rexact(r; T) = Qv^ + vi\ Pexact(q, s; T)|s=o , (26) 
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where specializing to the case of the 2D HO, 



Pcxact 



q,s;T) = -5^5^(-l)'^L„(2g2)Lfc(sV2)e- 



(g2+s2/4). 



1 



exp( "+y"^ ) + 1 



(27) 



n=0 fc=0 

is the ejcac? finite-temperature first-order density matrix. Putting s = in Eq. (27) gives the exact 
finite-temperature particle density for the 2D HO, viz., 



Pcxact 



oo oo ^ 



exp( "+y-^ ) + 1 



(28) 



n=0 fc=0 

In Table I below, we present a numerical comparison of the kinetic energies as obtained from 



-^exact 
-^TFT 
i^ETFT 



J re.act(r)ciV , c.f.Eq. (23) 

j TrFTiPexactlfi^r, c.f. Eq. (12) 
j TETFT[Pe^!,ct]d'^r, C.f. Eq. (12) 



(29) 
(30) 
(31) 



for N = A2 particles, with Table II providing the same calculation for = 420 particles. 



T 


-f^cxact 


-^TFT 


Ketft 


0.5 


93.8984 


93.4202 


93.6977 


0.8 


97.8489 


97.3152 


97.7439 


1.0 


101.3291 


100.7624 


101.2613 


2.0 


125.9995 


125.3372 


125.9854 


3.0 


157.9066 


157.2537 


157.9007 


4.0 


193.5600 


192.9616 


193.5591 


5.0 


231.2762 


230.6358 


231.1417 



Table I: A comparison of the total kinetic energy at various temperatures as determined from Eqs. (29-31) 
with iV = 42 particles. All quantities are measured in scaled units, as discussed in the text. The largest 
relative percentage errors in our tabulated data are AKtft ~ 0.5% and Ai^'sTFT ~ 0.2%. 



It is clear that the TFT KE, /Ctftj is always lower than the exact KE density, iCcxactj while the 
gradient corrections serve to improve the agreement with the exact result. This makes sense given 
that the gradient corrections take into account the curvature of the system imposed by the external 
potential, thereby increasing the kinetic energy of the system. It is nevertheless quite surprising 



T 


K 

-''-exact 


-^TFT 


Ketft 


0.5 


2879.2573 


2877.8822 


2878.3566 


0.8 


2892.3488 


2890.9642 


2891.7668 


1.0 


2904.3701 


2902.9579 


2903.9484 


2.0 


3002.5299 


3000.9349 


3002.4338 


3.0 


3158.2817 


3156.5131 


3158.2465 


4.0 


3361.6954 


3359.7806 


3361.6754 


5.0 


3603.0235 


3601.0015 


3603.0076 



Table II: As in Table I, but with N = 420 particles. The largest relative percentage errors in our tabulated 
data are AKtft ~ 0.06% and Ai^ETFT ~ 0.03%. 

how well the TFT functional does in describing the kinetic energy of the strongly inhomogeneous 
2D HO at finite temperature, even for small particle numbers. From our numerical calculations, we 
observe that a ten-fold increase in the number of particles reduces the largest relative percentage 
error by roughly a factor of ten; the better agreement between the TFT, ETFT and the exact KE, is 
in keeping with the expected result that in the large- limit, the TF approximation becomes exact. 

In Figs. [Hand [2l we present the KE densities, 7^xact(r), 7tft[p], and 7etft[p] with N = A2 
particles, at T = 0.2 and T = 2, respectively. As in Tables I and II, the exact spatial density, 
Eq. (28), has been used as input for the KE functionals. We have focused on a small number of 
particles, viz., N = 42, since deviations between the exact, TFT, and ETFT densities are more 
pronounced for N < O(IO^), particularly at low temperatures. 

Figure [H reveals several interesting aspects of the level of approximation at low temperatures. 
First, we note that the TFT (dot-dashed curve, black online) and the exact KE density (solid curve, 
red online) are almost indistinguishable within < r < 3. However, near the tail region (see figure 
inset), it is clear that the two KE densities are quite different; the TFT density is strictly positive 
definite, whereas the exact KE density falls below zero, before coalescing with the TFT density for 
r > 4.5. The ETFT KE density (dashed curve, blue online), on the other hand, does a relatively 
poor job of quantitatively capturing the behaviour of the exact KE density in the interior, but for the 
tail region, more accurately reproduces the exact result. Indeed, for r > 3.5, the ETFT and exact 
KE densities are indistinguishable on the scale of the inset. Moreover, in spite of the differences 
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Figure 1: (colour online) A plot of the exact, Eq. (26), TFT, i.e., the first term in Eq. (16), and the ETFT, 
Eq. (16), KE densities for N = A2 particles and T = 0.2. The solid (red online) curve is the exact KE 
density, the dot-dashed (black online) curve is the TFT KE density, and the dashed (blue online) curve is 
the ETFT KE density. Inset: Magnification of the tail region, where deviations between the three curves are 
most pronounced. The arrow indicates the classical turning point. Scaled units are used, as discussed in the 
text. 

between the exact and ETFT densities for r < 3, the ETFT is still a better approximation for the 
total KE, as evidenced by the data presented in Table 1. Therefore, while the gradient corrections 
are important for improving the total (i.e., integrated) KE, they are essential for describing the 
correct low-temperature behaviour of the exact KE past the classical turning point. 

Figure [21 presents the same data as in Fig. [U but with T = 2.0. At this temperature, the 
shell oscillations are already completely washed out, and the exact and ETFT KE densities are 
indistinguishable from each other, including the tail region; in the inset, there are actually three 
curves plotted, but the difference between the solid (red online) and dashed (blue online) curves 
cannot be resolved. Thus, at temperatures for which the shell effects are absent (i.e., T > 1), 
the ETFT is an excellent approximation to the exact finite-temperature KE density. Near the tail 
region, we see that the TFT KE density (dot-dashed curve, black online) is consistently too large, 
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r 



Figure 2: (colour online) As in Fig.[T] but with T = 2.0. Note that at this temperature, the shell oscillations 
are already completely washed away. Inset: Magnification of the tail region, clearly illustrating that the 
exact (solid curve, red online) and ETFT (dashed curve, blue online) curves are indistinguishable on the 
scale of the plot. The arrow indicates the classical turning point. 

thereby emphasizing the importance of the gradient corrections for a faithful description of the 
local behaviour of the exact KE density, even for small particle numbers. 

Finally, in Fig. |3l we illustrate the temperature and spatial dependence of the vW coefficient, 
/2(r;T) given by Eq. (18), for = 420 particles, and N = 42 (inset). As the temperature is 
increased, we see that the vW coefficient approaches the constant value /2(r; T) = 1, confirming 
our analytical results above for the Boltzmann regime. Figure [3] establishes that T > 18 is a 
sufficiently high temperature for the = 420 particle system to be treated as Boltzmann gas. The 
inset illustrates the expected result that for smaller particle numbers, one enters the Boltzmann 
regime at much lower temperatures. 
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Figure 3: A plot of the coefficient of the von Weizsacker term, Eq. (18), for N = 420 particles. The curves 
correspond, from lowest to highest, to T = 2 and T = 18, respectively, in steps of AT = 2. Note that by 
T = 18, /2(r; T) is already approaching its Boltzmann value, /2(r; T) = 1, represented by the dashed line. 
Inset: As in the main figure, but for iV = 42 particles at temperatures T = 1, 3, 5, 7. Scaled units have been 
used, as discussed in the text. 

IV. CONCLUSIONS AND OUTLOOK 

We have provided a formal justification for the inclusion of gradient corrections to the 2D 
KE density functional of an ideal Fermi gas at finite temperatures. Our numerical calculations 
have examined the quality of the TFT and ETFT functionals by comparing them against exact, 
analytical results for the 2D HO potential. We find that gradient corrections lead to an improved 
agreement for total KE when compared to the TFT approximation, and are necessary to correctly 
reproduce the quantum mechanical tunneling into the classically forbidden region exhibited by the 
exact KE density. Unfortunately, the non-analytic behaviour of the T = semiclassical densities 
at the classical turning point cannot be remedied within the present formalism. 

An extension of this work would be to develop the finite-temperature Dirac exchange func- 
tional, which could then be used in a fully self-consistent, finite-temperature Thomas-Fermi-Dirac 
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von Weizsacker (TFDW) DFT calculation similar to what has already been done at T = for 
low-dimensional Fermi systems.— It would be interesting to see whether the optimal, ad hoc, 
T = vW coefficient of 1 /8 could be motivated from a finite-temperature self-consistent TFDW 
calculation. 

Finally, we wish to point out that the results presented here may also find relevance in cur- 
rent experiments on ultra-cold, trapped Fermi gases, in which inter-atomic interactions may be 
tuned from essentially zero, to very strong, via the Feshbach resonance.— It is possible that the 
low-temperature shell oscillations, and their suppression as the temperature is increased, may be 
directly observable in cold atoms experiments on low-dimensional systems. 
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